QCD matter within a quasi-particle model and the critical end point 
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We compare our quasi-particle model with recent lattice QCD results for the equation 
of state at finite temperature and baryo-chemical potential. The inclusion of the QCD 
critical end point into models is discussed. We propose a family of equations of state 
to be employed in hydrodynamical calculations of particle spectra at RHIC energies and 
compare with the differential azimuthal anisotropy of strange and charm hadrons. 

1. QUASI-PARTICLE MODEL 

The equation of state (EoS) of strongly interacting matter is a central issue to un- 
derstand the hydrodynamical evolution in high-energy heavy-ion collisions. Ab initio 
evaluations based on QCD can deliver the EoS, e.g., in the form [Q] 

P = T 4 jco(T) + c 2 (T) (|) 2 + c 4 (T) (f) 4 + c 6 (T) ^ + ••■}. (1) 

To arrive at a flexible parametrization of the EoS we have developed a quasi-particle model 
[121 El using dynamically generated self-energies and a non-perturbative effective coupling 
as essential input. A comparison of the model with recent lattice data is exhibited in Fig. 1. 
The peak of C4 and the dipole shape of Cq emerge from a change of the effective coupling 
at the pseudo-critical temperature T c (/i B = 0) dictated by the data. The calculation of 
the coefficients 02,4,6 requires the knowledge of the pressure p(T,fi B ) at finite values of 
baryo-chemical potential which is delivered by Peshier's equation [|2J. The baryon 
number susceptibility follows directly from p(T,fi B ), see Fig. 1. The peak developing for 
\Ib > 300 MeV can be considered as a signal of some critical behavior. 

2. CRITICAL END POINT 

QCD for finite quark masses displays a critical end point (CEP) [ 0J where a phase 
transition of first order sets in at finite values of \xb- The CEP belongs to the universality 
class of the 3D Ising model characterized by a set of critical exponents. A convenient 
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Figure 1. The Taylor expansion coefficients Co 5 2,4 (top row, horizontal lines depict the 
corresponding terms in the bag model) and cq (left panel in bottom row) of the equation 
of state as a function of the scaled temperature. QCD lattice data from [HJE]. Right 
panel in bottom row: baryon number susceptibility (//# = 450, 330, 300, 150 MeV from 
top to bottom). For details see [E]. 



parametric representation of the singular part (formulated in terms of the entropy density 
s) of the EoS can be found in [E]. The complete EoS is accordingly 

S Sreg ^sing- (2) 

While the location of the CEP is determined by lattice QCD calculations [|7j the size of 
the critical region is barely known. A simple toy model demonstrating the effect of s S i ng 
for given s reg is s reg = 4coT 3 + 2c2^ 2 B T and s S i ng = s reg A tanh5 c with the same S c as 
in [E| and c = (32 + 21iV / )7r 2 /180 and c 2 = N f /18 for N f = 2. We assume a critical 
line given by T(fi B ) = T c (l - c(/i B /T c ) 2 ) with c = 0.0077 from [ 9J and locate the CEP 
at //£ = 330 MeV according to [Ej. This information is needed to map the coordinates 
of s sing into the T — ji B plane. A typical pattern of isentropic trajectories is exhibited 
in Fig. 2. The pattern (left panel) seems to be generic (see also [|H] for another model 
with CEP): trajectories which are left [right] to the CEP for small values of the strength 
parameter A are attracted [repelled], but at variance to lattice QCD results (right panel). 

The proper implementation of CEP phenomena in our quasi-particle model is under 
consideration. 
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Figure 2. Isentropic trajectories and the (pseudo) critical line in the T — fi B plane. Left 
panel: a toy model including the CEP (entropy per baryon = 100, various strength 
parameters A as indicated). Right panel: quasi-particle model adjusted to lattice QCD 
data [ HI |5] via c (solid curves) or c 2 (dashed curves); both fits agree within 10% for 
p(T,fi B = 0); entropy per baryon = 200, 100, 50, 33, from left to right. 



3. A FAMILY OF EQUATIONS OF STATE, HYDRODYNAMICS AND 
AZIMUTHAL ANISOTROPY 

It is expected that lattice QCD data at high temperature are realistic, as improved 
actions have been used. At low temperature, the employed quark masses are too large 
to give realistic results. These however agree with the resonance gas model once analog 
assumptions are implemented [HD|- It is therefore reasonable to use the resonance gas 
model at low T and employ our quasi-particle model to extrapolate lattice QCD results 
both to larger fis and smaller T as long as a systematic chiral extrapolation can not be 
done safely due to lacking lattice QCD input. 
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Figure 3. A family of EoS's by combining our quasi-particle model and the resonance gas 
model for baryon densities relevant for top RHIC energies: pressure (left, QPM(4.0), 
QPM(2.0), QPM(1.25), QPM(l.O), from top to bottom) and squared sound velocity 
(right) as a function of energy density. Bag model results are depicted by dashed curves. 
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We generate a family of EoS's by keeping the matching point to the resonance gas EoS 
fixed and interpolate linearly to a given matching point (given by energy density which 
serves as label) of the quasi-particle model. Such a family is exhibited in Fig. El Surprising 
is the similarity of QPM(l.O) and the bag model EoS used in [EJ. Equipped with such a 
QCD based EoS one can compare results of hydrodynamical calculations with data, e.g., 
for the azimuthal anisotropy v-i. Examples are exhibited in Fig. 0] for A, H, Q, <f), D. 




Figure 4. v-i as a function of transverse momentum. Left: strange baryons (solid curves 
QPM(l.O), QPM(2.0), QPM(4.0) [from top to bottom, with decoupling temperatures 
adjusted to data up to 2.5 GeV/c] and squares: A, dashed curve [QPM(2.0] and stars: S, 
dotted curve [QPM(2.0)] and circles: Q, data from [Q2]), right top: <fi (data from [IT5]). 
right bottom: D (data from [|14|). both for QPM(2.0). Impact parameter b = 5.2 fm. 
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